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Abstract 

The influence of excited-state absorption (ESA) and two-exciton processes on a co- 
herent population transfer with intense ultrashort chirped pulses in molecular systems 
in solution has been studied. An unified treatment of adiabatic rapid passage (ARP) 
in such systems has been developed using a three-state electronic system with relax- 
ation treated as a diffusion on electronic potential energy surfaces. We have shown 
that ESA has a profound effect on coherent population transfer in large molecules that 
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necessitates a more accurate interpretation of experimental data. A simple and phys- 
ically clear model for ARP in molecules with three electronic states in solution has 
been developed by extending the Landau-Zener calculations putting in a third level 
to random crossing of levels. A method for quantum control of two-exciton states in 
molecular complexes has been proposed. 

1 Introduction. 

The possibility of the optical control of molecular dynamics using properly tailored pulses 
has been the subject of intensive studies in the last few years [U(2l[3llll[5l[6],[71[8|,[9l EH [TTJ, 
H21 H31 HH H5l HH HZl HH HH SOI [2U [22I Eg. Chirped pulses can selectively excite coherent 
wave packet motion either on the ground electronic potential energy surface of a molecule or 
on the excited electronic potential energy surface due to the intrapulse pump-dump process 
ID El ED E2]- In addition, they are very efficient for achieving optical population transfer 
between molecular electronic states. Total electronic population inversion can be achieved 
using coherent light-matter interactions like adiabatic rapid passage (ARP) in a two- or three- 
state system [2U [22], which is based on sweeping the pulse frequency through a resonance. 

Since the overwhelming majority of chemical reactions are carried out in liquid solution, 
adiabatic passage in molecules in solution was studied for two-state electronic system (ARP) 
in Refs.[2l)l 1271 [28] . and for stimulated Raman adiabatic passage (STIRAP) configuration 
in Refs.|29[ [30] . It has been shown in Ref.[26] that relaxation does not hinder a coherent 
population transfer for positive chirped pulses and moderate detuning of the central pulse 
frequency with respect to the frequency of Franck-Condon transition. 

However, a two electronic state model for molecular systems is of limited utility. Indeed, 
excited-state absorption (ESA) occurs for majority of complex organic molecules [3TI [25]. 
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Even a molecular dimer consisting from two-level chromophores has an additional excited 
state corresponding to two-exciton excitation. An unified treatment of ARP in such sys- 
tems can be developed using three-state electronic system interacting with reservoir (the 
vibrational subsystems of a molecule (chromophores) and a solvent). 

More often than not ESA in complex organic molecules corresponds to a transition from 
the first excited singlet state Si to a higher singlet state S n (n > 1), which relaxes back 
to Si very fast [32], [331 EH ED]. Therefore, it would look as if ESA does not influence on 
population transfer So — > Si from the ground state So. However, in the presence of ESA an 
exciting pulse interacts with both So — > Si and Si — > S n transitions. It is well known that 
coherent optical interactions occurring in adjacent optical transitions in a three-state system 
markedly affect each other. The examples are STIRAP, lasing without inversion, coherent 
trapping, electromagnetically induced transparency and others. (For textbook treatments of 
these effects see, for example, |35j). Therefore, one would expect an appreciable change of 
a population transfer S — > Si with chirped pulses in the presence of excited state absorp- 
tion in the coherent regime when the chirp rate in the frequency domain is not large and, 
consequently, the pulse is rather short. 

Our objective is to answer the following questions: "How do ESA and two-exciton pro- 
cesses influence on a coherent population transfer in molecular systems in solution? What 
is the potential of chirped pulses for selective excitation of the single and two-exciton states 
and their selective spectroscopy?" 

In addition, the three-state system under discussion enables us to consider STIRAP as 
well. Therefore, we shall also briefly concern slowing down the pure dephasing on STIRAP 
in intense fields when relaxation is non-Markovian. 

The outline of the paper is as follows. In SecfSJwe present equations for the density matrix 
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of a three-state molecular system under the action of shaped pulses when the interaction 
with a dissipative environment can be described as the Gaussian-Markovian modulation (so 
called the total model). In Secj3]we formulate a number of approaches to this model that 
enables us, first, to clarify the underlying physics and, second, to understand the validity of 
the results obtained by the total model. The ESA effects on ARP in complex molecules are 
considered in SecJU In Secj5] we study population transfer in molecular dimers with taking 
into account two-exciton processes. In Secj6]we consider slowing down the pure dephasing 
on STIRAP in strong fields when the system-bath interaction is not weak (non-Markovian 
relaxation). We summarize our results in SecJTJ In the Appendix we extend calculations 
of two-photon excitation of a quantum ladder system by a chirped pulse [36J to non-zero 
two-photon detuning. 



2 Basic equations 

Let us consider a molecular system with three electronic states n = 1,2 and 3 in a solvent 
described by the Hamiltonian 
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H = \n) [E n + W n (Q)] (n\ (1) 

n=l 

where E 3 > E 2 > Ei, E n is the energy of state n, W n (Q) is the adiabatic Hamiltonian of 
reservoir R (the vibrational subsystems of a molecular system and a solvent interacting with 
the three-level electron system under consideration in state n). 

The molecular system is affected by two shaped pulses of carrier frequencies u\ and uo 2 



E(t) = - E i(t) + c - c - = 2 Yl ^ (*) 6X pHW + m (*)] + c - c - ( 2 ) 

i=l,2 i=l,2 

which are resonant to optical transitions 1^2 and 2^3, respectively (ladder configura- 



tion). Here Si (i) and ^ (t) describe the change of the pulse amplitude and phase, respec- 
tively, in a time t. The instantaneous pulse frequencies are Ui (t) = u)i — . 

The influence of the vibrational subsystems of a solute and a solvent on the electronic 
transition can be described as a modulation of this transition by low frequency (LF) vibra- 
tions {u s } [37J 38J. In accordance with the Franck-Condon principle, an electronic tran- 
sition takes place at a fixed nuclear configuration. Therefore, for example, the quantity 
w(Q) = W 2 (Q) - Wi(Q) - (W 2 (Q) - Wi(Q))i is the disturbance of nuclear motion under 
electronic transition 1 — > 2. Here () n = Tr R (...pji n ) denotes the trace operation over the 
reservoir variables in the electronic state n, 

p Rn = exp {-pW n ) /Tr R exp (~pW n ) , (3 = l/k B T. 

The relaxation of electronic transition 1—^2 stimulated by LF vibrations is described by 
the correlation function K(t) = (u(0)u(t)) of the corresponding vibrational disturbance with 
characteristic attenuation time r s [12], [38] - We suppose that ftw s k R T. Thus {oa,} is an al- 
most classical system and operators W n are assumed to be stochastic functions of time in the 
Heisenberg representation. The quantity u can be considered as a stochastic Gaussian vari- 
able. We consider the Gaussian-Markovian process when K(t)/K(0) = S(t) = exp(— |i|/r s ). 
The corresponding Fokker-Planck operator Lj = r~ 
diffusion in the effective parabolic potential 

^ (q) = E 3 + l -u 2 (q - d,f (3) 

of electronic state j where r" 1 = D n j3Co 2 and D is the diffusion coefficient. Going to a di- 
mensionless generalized coordinate x = quj\fj3, one can obtain the equations for the elements 
of the density matrix pij(x,t) by the generalization of the equations of Ref. [26J. Switching 
to the system that rotates with instantaneous frequency 



i a 2 

f3u> 2 dq 2 



+ (9-<*i)£ + l 



describes the 
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p 12 (x,t) = pi 2 (x,t) exp[-i(cjit - Vi (*))], P2z{x,t) = p 23 (x,t) exp[-i(u 2 t - V?2 (*))], 
Pi 3 OM) = P13OM) exp{-i[(cji + cj 2 )t - + V?2 (*))], (4) 



we get 



<9 

— pn(ar, i) = Im[fiipi 2 (x, *)] + LiPn(a;, t) 
at 

<9 

— p 22 (a:, = - Im[fiipi 2 (x, t) + fi 2 p 32 (x, *)] + L 2 p 22 (x, t) + 2r 32 p 33 (x, f) 
or 
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— p 33 (x, t) = - Im[fi 2 p 23 (:r, t)\ + (L 3 - 2r 32 )p 33 (:r, i) (5) 
— pi 2 (x,t) = i [u 21 - ui(t) - (hP)~ 1 x 2 x\ pn(x,t) + -Ql[p 22 (x,t) - p u (x,t)]- 

- ^2Pis(x, t) + L 12 p l2 (x, t) (6) 

d % 

Qj.M x > *) = * [ w 3i - - wi(t) - (h(3)~ l x 3 x\ p 13 (x, t) + -ttlp 23 (x, t)- 

1 

- ^* 2 Pi2(x, t) + (Li 3 - r 32 )pi 3 (:r, i) (7) 
g^P^( x , t) = i [(^3i - w 2 i) - w 2 (t) - (a* - a: 2 ) x] p 23 (:r, i) + -fi 2 (t) [p 33 (x, t) - p 22 (:r, t)] + 

+ ^ipi 3 (x, t) + (L 23 - r 32 )p 23 (a;, i) (8) 

where Vt\ = D 2 \E\jfi and Q 2 = D 32 £ 2 /h are the Rabi frequencies for transitions 1 — > 2 
and 2 — >• 3, respectively. Here u^i = cofl + xf/(2h(3) is the frequency of Franck-Condon 
transition 1 — > i, cjfj = (Ej — Ej)/h is the frequency of purely electronic transition j — > i, 
Djj are matrix elements of the dipole moment operator, 2r 32 is a probability of nonradiative 
transition 3 — > 2 for the excited state absorption problem (see below); \xj\ = {h^u^J.) 1 / 2 is a 
dimensionless shift between the potential surfaces of states 1 and j (xi = 0), which is related 
to the corresponding Stokes shift u)\{ of the equilibrium absorption and luminescence spectra 
for transition 1 — > j. The last magnitude can be written as u) l s { = h/3aH where a^ s denotes 



the LF vibration contribution to a second central moment of an absorption spectrum for 
transition 1 — > j. The terms 

i, = v 1 (^ + («-*,)s + 1 ) < 9 > 

on the right-hand side of Eqs.fjS]) describe the diffusion in the corresponding effective parabolic 
potential 

U j (x) = E j + ±(x-x j ) 2 (j = 1,2,3), (10) 

Lij = [Li + Lj)/2. 

The partial density matrix of the system (x, t) describes the system distribution with a 
given value of x at time t. The complete density matrix averaged over the stochastic process 
which modulates the system energy levels, is obtained by integration of p^ (x, t) over the 
generalized coordinate x: 

(p)ij (t) = J Pij (x, t) dx (11) 

where diagonal quantities (p)jj (t) are nothing more nor less than the populations of the 
electronic states: (p)jj (t) = rij, ri\ + n 2 + n 3 = 1. 

We solve coupled Eqs.©-®, using a basis set expansion with eigenfunctions of diffusion 
operator L 13 , similar to Ref. [2"o] . 

The solutions, corresponding to the procedure described in this section, are termed the 
total model for short, bearing in mind that they take into account all the relaxations (diffu- 
sions) related to populations and electronic coherences between all the electronic states. 

3 Approximate models 

In this section we describe a number of approaches to the total model (Eqs.©-©). 
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3.1 System with frozen nuclear motion 

For pulses much shorter than r s one can ignore all the terms ~ Li , Lij on the right-hand 
sides of Eqs.(j5])-(J8]). It means that our system can be described as an ensemble of inde- 
pendent three-level systems with different transition frequencies corresponding to a pure 
inhomogeneously broadened electronic transitions. In this case the density matrix equations 
can be integrated independently for each x. After this the result must be averaged over x. 
Solutions of the undamped equations for the density matrix are interesting from the point 
of view of evaluation of the greatest possible population of excited states due to coherent 
effects, because these solutions ignore all the irreversible relaxations destructing coherence. 
In addition, a comparison between the latter solutions and calculations for the total model 
enables us to clarify the role of relaxation in the chirp dependence of population transfer 
(see SecJH below) . The approach under discussion in this section is termed "relaxation- free" 
model for short. 

3.2 Semiclassical (Lax) approximation 

For broad electronic transitions satisfying the "slow modulation" limit, we have c^t^ 3> 1, 
where cr^, is the LF vibration contribution to a second central moment of an absorption 
spectrum for transition % — > j. In the last case electronic dephasing is fast, and one can 
use a semiclassical (short time) approximation [39]. This limit is also known as the case of 
appreciable Stokes losses because the perturbation of the nuclear system under electronic 
excitation i — > j (a quantity Wj — Wi) is large. Then one can ignore the last term LijPij(x, t) 
on the right-hand side of the corresponding equation for the nondiagonal element of the 
density matrix [251 EH EE2J HI] that describes relaxation (diffusion) of pij(x,t) (Eqs.§6$) and 
([8])). The solutions, which correspond to disregarding terms Lijpij(x,t) for broad electronic 
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transitions i — > j are termed "partial relaxation" model for short It is worthy to note 
that the "partial relaxation" model offers a particular advantage over the total model. The 
point is that the first can be derived not assuming the standard adiabatic elimination of 
the momentum p for the non-diagonal density matrix [JT], which is incorrect in the "slow 
modulation" limit [12]. This issue is quite important in the light of the limits imposed on 
Eqs.((6]) and ([8]) for nondiagonal elements of the density matrix [431144]. 

Indeed, in the Wigner representation (15J HfJ H7] equation for p~i 2 may be written in the 
rotating frame as (see Eq.(j6])) 



i 

■^Pwu{q,P,t) = i[(U 2 (q) - Ui (q))/h-aji(t)}p wi2 (q,p,t) - -fi 2 p wi3 (g,p, t) + 

+ ^*[PW22(q,P,t) - Pwn{q,P,t)] + L FP12 p wl2 (q,p,t) (12) 

Eq. f[T2l has been derived for harmonic potentials, Eq.([3]), by generalization of equations of 
Refs. gg EH S21 E] where 



d d 
Lfpi2 = -Ptt + 

oq op 



is the Fokker-Planck operator for overdamped Brownian oscillator with attenuation constant 
7- 

In the case of appreciable Stokes losses when the perturbation of the nuclear system 
under electronic excitation 1 — > 2 (a quantity (U 2 (q) — U\ (q))/h — u> 2 [) is large, the quantity 
Pwi2{q,P,t) oscillates fast due to the first term on the right-hand side of Eq.(fl2l) (see also 
Ref . [42] ) . Therefore, to the first approximation, on can neglect changes of Pwi2(q,p,t) due 
to the last term on the right-hand side of Eq. fll2l) . Neglecting this term, integrating both 
side of Eq. fTl~2|) over momentum, and bearing in mind that 
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Pij(q,t) = / Pwij(q,P,t)dp 



(13) 



and x = qujy/]3, we get 



— pi 2 (x,t) = ?[^2i-^i(i)-(^/9) _1 ^2x]pi2(s,t) + -finP22(a;,i)-pii(a;,t)]--fi2Pi3(a;,t) (14) 

that is nothing more nor less Eq.([6]) without the last term Li2p~i2(x,t) on the right-hand 
side. As a matter of fact, a derivation of Eq. ffT4"]) does not involve the assumption that the 
momentum is instantly equilibrated. The same can be done with Eq.® for p 2 3- 



4 Adiabatic population transfer in the presence of excited- 
state absorption 

We shall study the ESA effects on ARP in complex molecules by the example of Coumarin 
153 in liquid solution [31J. In the frequency domain, the electric field can be written as 
|-E(£>) | exp[i<&(a))] and the phase term $(£>) can be expanded in a Taylor series <&(£?) = 
$(o;) + (1/2)$" (u) (£u — uj) 2 + ... We shall consider linear chirped pulses of the form 

E(t) = S exp[-i(5 2 - i»(t - t ) 2 ] (15) 
where the parameters 8 and p are determined by the formulae [TT1 [12] : 

5 2 = 2{r 2 + [2$" (co) /r p0 ] 2 }~\ I* = -4*" M K + 4$" 2 («)] ^ , (16) 

Tpo = tpo/V2 In 2, t p o is the pulse duration of the corresponding transform-limited pulse. Fig{T] 
shows populations of electronic states after the completion of the one pulse action as functions 
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-5 -4 -3 -2 -1 1 2 3 4 5 
<D"(v) (x10 4 fS 2 ) 

Figure 1: Populations of electronic states after the completion of the pulse action as functions 
of $"(z/) in a three-state system. Calculations without decay of the upper state 3 into state 
2: rii (dotted line), n 2 (solid line), (dashed line). Line with hollow circles - n 2 in the 
model with fast decay 3^2 r 32 =10 ps _1 . For comparison we also show n 2 for a two-state 
system (line with squares). Total relaxation model with diffusion of all matrix elements. 

of the chirp rate in the frequency domain <§>"{y) = 47r 2 $" (uj). For the molecule under 
consideration a two-photon resonance occurs at the doubled frequency of the Franck-Condon 
transition 1 — ► 2. Absorption spectrum corresponding to transition 1 — ^ 3 is rather narrow 
that means x 3 = 0. The values of parameters for Fig{T]were as follows: the pulse duration 
of the transform-limited (non-chirped) pulse t p0 = 10 fs, uj\ 2 = 2686 cm" 1 , D u = D i2 = 6 
D [31], t s = 70 fs, the saturation parameter, which is proportional to the pulse energy 
[26] . Q' = y/7r\Di 2 S max \ 2 t p / (2h 2 a/2o^J) = 5; the one-photon resonance for Franck-Condon 
transition 1 — > 2 occurs at the pulse maximum, i.e. uj = u 2 i- 

FigfS] contrasts calculations using the total model (FigJTJ) with those of the partial relax- 
ation model. The latter includes both diffusion of all the diagonal elements of the density 
matrix and one off-diagonal element pi%. The point is that transition 1 — > 3 occurs with- 
out changing the state of vibrational subsystems of a molecule and a solvent, and therefore 
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-5 -4 -3 -2 -1 1 2 3 4 5 
<D"(v) (x10 4 fS 2 ) 



Figure 2: Populations of electronic states n\ (dotted lines), n 2 (solid lines) and n 3 (dashed 
lines) after the completion of the pulse action calculated without decay of the upper state 3 
into state 2 as functions of <&"{y). The partial relaxation and the total models - lines with 
and without hollow circles, respectively. All the parameters are identical to those of Fig{TJ 



can not be described in a semiclassical (short time) approximation. Figj2] shows a good 
agreement between calculation results for the models under consideration. 

One can see from Fig{TJ first, that population 712 for a molecule with a fast decay 3^2, 
which closely resembles experimental data [TTJ for LD69Qj, is distinctly different from that 
of a two-state system for |$"(z/)| < 15 • 10 3 fs 2 when the excited pulse is rather short. This 
means that the excited state absorption has a profound effect on coherent population transfer 
in complex molecules. Second, n 3 strongly decreases when |$"(z/)| increases. 

To understand these results, we will consider first two transitions separately. One can 
obtain the following criterion for the adiabaticity of one transition in the absence of relax- 
ation: Q' » 1 where Q' is the saturation parameter. It conforms to the value of Q 1 = 5 
used in our calculations. The condition Q' » 1 follows from the adiabatic criterion for a 
two-level system: 



According to Ref . |23j , LD690 shows ESA. 
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Figure 3: Populations of electronic states n\ (dotted line), n 2 (solid line), n 3 (dashed line) 
and n 2 + n 3 (line with hollow circles) after the completion of the pulse action as functions 
of <&"(u) for the relaxation-free model t s — > oo. Other parameters are identical to those of 
Fig{TJ In the case under consideration the combined population n 2 + does not depend on 



^|«i^wr m 

where fi 12 (t) = \D 2 i^ 32 S(t)\/h are the Rabi frequencies for transitions 1 — > 2 and 2^3, 
respectively. Adiabatic criterion Eq. (fT7|) was fulfilled in our simulations for both transitions 
1 — > 2 and 2 — * 3 at any <&"(y). However, Fig{T] shows that strongly decreases when 
|$"(z/)| increases. To clarify the reasons for strong decreasing n 3 it is instructive to carry 
out the corresponding calculations for the relaxation-free model of Sec J3. II shown in Figj3l 
In this case excitation of state 3 with a transform-limited pulse is slightly more effective as 
compared to a strongly chirped pulse of the same energy. The point is that a two-photon 
resonance occurs for a number of spectral components of a transform-limited pulse and only 
at the maximum of a strongly chirped pulse. However, Figfjdoes not show strong decreasing 
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the population of state 3 when $"(z/) increases. This means that relaxation is responsible 
for strong decreasing n% as a function of $"(z/) in spite of the fact that relaxation does not 
destroy ARP when the Rabi frequencies exceed the reciprocal irreversible dephasing time 
(T') _1 |2T] 



fti, 2 » l/T' (18) 

The last condition was fulfilled in our simulations at least for |<&"(i/)| < 10 4 fs 2 . 

To clarify this issue, we shall consider a population transfer between randomly fluctuating 
levels. 



4.1 Population transfer between randomly fluctuating levels 

The picture of randomly fluctuating levels |27J offers a simple and physically clear explanation 
of numerical results [26J obtained for population transfer in a two-state system. Here we shall 
generalize the Landau-Zener (LZ) calculations putting in a third level [50J to random crossing 
of levels. 

Let us write the Schrodinger equations for the amplitudes of states 01,2,3 for the system 
under consideration. Switching to new variables a k '- 



a k = a k exp 



we obtain in the rotating wave approximation 
/ \ / 



- / U 2 dt 



(19) 



.d 
dt 



a 1 
a 2 



V 



(Ux-U 2 )/h + wi(t) -fii/2 

-fti/2 -0 2 /2 

-n 2 /2 {u 3 -u 2 )/h-u 2 (t) 



\ ( - \ 



a 2 



\ a3 J 



(20) 



Throughout this section effective parabolic potentials ffTUl) are considered as functions of gen- 
eralized coordinate a = x^faf s -uj X2 : Uj(a) = Ej+^ia+^u^l^u^+^-l) 39 ^ \[u^ t }} 2 '. 
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Here 

(U 3 - U 2 )/h - u 2 {t) = [{uof 2 + u%/2) - uo 2 ] + o + fi 2 t (21) 
for X3 = X\ (that corresponds to Coumarin 153), and 

(U x - U 2 )/h + Ul (t) = [u> x - « - lu^/2)] + o - pit, (22) 

for linear chirped pulses u 12 (t) = Ui <2 — \i\ i2 t. 

Let us define instantaneous crossings of state 2 with photonic repetitions 1' and 3' of 
states 1 and 3, respectively. They are determined by the conditions that quantities Eqs.(!2~TT) 
and ( 122]) are equal to zero: 



ai 2 (t) = Ki - u%/2) - wi + Hit = o 12 (0) + Hit (23) 
a 23 (t) =U2- (^32 + /2) - Vit = o 23 (0) - Hit 

Near the intersection points one can consider a as a linear function of time. For small t, 
a(t) ~ 0:12(0) + at. Let 012(0) = 023(0), i.e. states 2, 1' and 3' cross at the same point when 
t = 0. This means 



U)f x + Ujf 2 = LOl + U) 2 



(24) 



i.e. the two-photon resonance occurs for t = 0. Then Eqs. ff20|) take the following form 



dt 



a 2 



ya 3 j 



\ 



{a-Hi)t 

-fix/2 -n 2 /2 

-Q 2 /2 (a + H 2 )t 



\ ( „ \ 



) 



a 2 



(25) 



that can be reduced to Eqs.(2) of Ref.|50j. Using the solution obtained in |50j and consid- 
ering identical chirps when Hi — ^2 = Hi we f° r the initial condition |oi(— oo)\ 2 = 1, 
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^2,31-00 J 







\a 3 (oo J 



(1-P)(1-Q) for -\fji\<a< \fji\ 
P(l — P)(l — Q) for both a > —fi when fi < 0, and d < — /x when /i > 
Q(l — P)(l — Q) for both a < fi when /x < 0, and a > fi when // > 



(26) 



where 



P = exp 



, Q = exp 



(27) 



4|d — n\J V 4|ck + /i| 

Similar to Ref.[2Z], we consider a as a stochastic Gaussian variable. Consequently, we 



must average Eqs. ff26|) over random crossing of levels described by Gaussian random noise 
induced by intra- and intermolecular fluctuations. It can be easily done for a differentiable 
(non-Markovian) Gaussian process [27], bearing in mind an independence of a and a from 
each other for such processes. Therefore, we shall consider in this section a differentiable 
( non-Markovian) Gaussian noise, as opposed to previuous sections. In addition, we consider 
a slow modulation limit when cx^t^ >> 1. Averaging Eqs. fl26p . we obtain the following 
expression for the population of state 3 when \i > 



n 3 = / da[ P{l-P)(l-Q)f{a,a)da+ (1 - P)(l - Q)f(a, a)da 



oo J — oo 
oo 



I Ml 



/■oo 

/ Q{l-P){l-Q)f{oc,a)da] 
J\iA 



(28) 



Here f(a, a) is the joint probability density for a and its derivative a: 



f(a,a) 



exp 



a 2 a 2 



+ 



2(J 2s 2fc(0) 



(29) 
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k (0) is the second derivative of the correlation function k(t) =< a(0)a(t) >= a^s ex P(~\t\/ T s) 
of the energetic fluctuations evaluated at zero. Eq. (l28|) is written for /i > (negatively 
chirped pulse). One can easily show that is symmetrical with respect to the chirp sign. 
The point is that a simple stochastic model of this Section misses any chromophore's effects 
on bath, in particular the dynamical Stokes shift (see Ref.[5T] for details). This is opposite 
to the models of previous sections, which do describe the dynamical Stokes by the drift term 
(the second term on the right-hand side of Eq.Q). 

Integrating Eq. (l28]) with respect to a and entering a dimensionless variable y = a / \fi\ , 
we get 



n 3 = {^\J_ 1 P (! " P )(! " Q) ex P ("f V 2 ) d V + / V " P )(! " Q)] ex P (-\ y2 ) dy 
+ Q(l - P)(l - Q) exp (-^ 2 ) dy] (30) 



where 



exp I ~%W=T\) and Q = exp ("25TH ' ' (31) 



and 



x= 4^ >0 , e = -^>0 (32) 
2|/i| fc(0) 1 J 

are dimensionless parameters. 

When adiabatic criterion Eq.( fl7l) is satisfied, parameter x is much larger than 1 since 
\duj(t)/dt\ = |/i| for a linear chirped pulse. Then the integrals on the right-hand side of 
Eq. (l30l can be evaluated by the method of Laplace, similar to Ref.[27j. The result is espe- 
cially simple for strong interaction, Eq. (I18j) . where the irreversible dephasing time of transi- 
tions 1 — > 2 and 2 — > 3 is given by [27] T' = l/[— A;(0)] 1//4 . Then, as one can see also from 
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Eqs.f l30p and (|31[) . the main contribution to n 3 is given by 

n ^{i L exp (~H dy = erf (^) (33) 

Since erf(1.5) = 0.966, we obtain that relaxation does not hinder a population transfer to 
state 3 when 



\n\T' 2 > 2 (34) 

For strongly chirped pulses [52], fi\T' 2 /V2 « 2 v / 2vr 2 T' 2 /|$"(z/)|. 

Eq. (134j) expresses an extra criterion for coherent population transfer to those we have 
obtained before for a two-level system [27], Eqs. ffTTl) and ffT8l) . New criterion (1341 implies 
conservation of the "counter-movement" of the "photonic repetitions" of states 1 and 3, 
in spite of random crossing of levels. Condition ( 134"1) is exemplified by FiglU In addition, 
FigJH shows an excellent agreement of simple formula (1331 with numerical calculations. It 
is worthy to note that condition (fl8|) was fulfilled in our simulations, though in the last 
case T' = (T s /a 2s ) 1 ^ 3 is determined independently of k(0) [S3], which does not exist for the 
Gaussian-Markovian process. 

4.2 Influence of excited-state absorption when detuning from two- 
photon resonance occurs 

For Coumarin 153 in liquid solution considered above a two-photon resonance occurs at 
the doubled frequency of the Franck-Condon transition 1 — > 2. In this section we consider 
populations of electronic states when the condition for two-photon resonance is violated. 
Figs. [5] and M show populations of electronic states for the total model after the completion 
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Figure 4: Population of state 3 as a function of the irreversible dephasing time T" for $"(z/) = 
10 4 fs 2 calculated by Eq. fl33l) (solid line with circles) and numerical solution of Eqs. ([HJ 
(dashed line with squares). n 3no relaxation = n 3(T' — > oo). Other parameters are identical to 
those of FigJH 




Figure 5: Populations of electronic states after the completion of the pulse action as functions 
of $"( v) in a three-state system. The frequency of purely electronic transition 3 — > 2, oj^i 
decreases by o;^ 2 /4 with the conservation of £3 = 0. Calculations without decay of the upper 
state 3 into state 2: rii (dotted line), n 2 (solid line), n 3 (dashed line). The corresponding 
populations in the model with fast decay 3^2 r 32 = 10 ps^ 1 are shown by the same lines 
with hollow circles. 
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Figure 6: Populations of electronic states after the completion of the pulse action as functions 
of $"(z/) in a three-state system. Equilibrium position of state 3 is offset to the right by 
^3 = ^2/2 and down so that frequencies of Franck-Condon transitions 1 — > 2 and 2 — > 3 are 
equal: U21 = u 32 - Calculations without decay of the upper state 3 into state 2: n\ (dotted 
line), n 2 (solid line), n 3 (dashed line). The corresponding populations in the model with fast 
decay 3^2 = 10 ps _1 are shown by the same lines with hollow circles. 
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of the pulse action as functions of $"(z/) for the same values of parameters as for Fig{T]with 
the only difference concerning the position of state 3. The frequency of purely electronic 
transition 3^2 uf 2 decreases by uj\^ /4 with the conservation of £3 = for Fig. [5j Equilib- 
rium position of state 3 is offset to the right by £3 = £2/2 and down so that frequencies of 
Franck-Condon transitions 1 — > 2 and 2^3 are equal: U21 = W32 for Fig. [6l 

One can see from FigsJIJ [5] and El first, that population n\ and, as a consequence, n 2 + n3 
depend only slightly on the occurrence of fast decay 3^2. Second, populations n 2 and 
ri3 in the absence of fast decay 3 — > 2 are very sensitive to the violation of the two-photon 
resonance condition. However, a behavior of n 2 , when fast decay 3 — > 2 occurs, and n\ 
as functions of $"(z/) is very similar for the figures under discussion, regadless of the two- 
photon resonance condition. Experimental measurements commonly correspond to n 2 and 
are carried out under the fast decay 3^2 conditions. Thus the behavior of n 2 for fast decay 
3^2 shown in FigsfjjJ [5] and [6] is rather versatile. 

5 Population transfer in the presence of two-exciton 
processes. Selective excitation of single and two- 
exciton states with chirped pulses 

Consider a dimer of chromophores each with two electronic states described by the Frenkel 
exciton Hamiltonian [5H [55l [56] and excited with electromagnetic field Eq.(j2J). The Hamil- 
tonian of the dimer is given by 

H= ^mB+B m + hJ{B+B 2 + B+B 1 )+H bath + H eb - D m • E(t)(5+ + fi m ) (35) 

771=1,2 771=1,2 
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where £> + = |ra)(0| (B m = |0}(m|) are exciton creation (annihilation) operators associated 
with the chromophore m, which satisfy the commutation rules [B n , £+] = 5 nm {l — 2B^ n B m ), 
5 nm is the Kroenecker delta; |0) and \m) denote the ground state and a state corresponding 
to the excitation of chromophore m, respectively. D m is the transition dipole moment of 
molecule m, Hb at h represents a bath and H e b its coupling with the exciton system. We 
assume that the bath is harmonic and that the coupling is linear in the nuclear coordinates 

Heb = -ft ^ a mnB+B n (36) 

mn 

where a mn represent collective bath coordinates. Ml^M^) and ft J are the exciton energy 
of 1 (2) chromophore and their coupling energy at the equilibrium nuclear coordinate of the 
ground electronic state. One can consider a mn as diagonal: a mn = a m 5 nm on the assumption 
that the electronic coupling constant fluctuation amplitude is negligibly smaller than the site 
energy fluctuation amplitude [35] . 

Diagonalizing the electronic Hamiltonian 



H e = h & m B+B m + hJ{B+B 2 + B+B 1 ) (37) 

m=l,2 

by unitary transformation [37] 



U 



-i 



/ \ 

cos 6 sin 6 



— sin 9 cos 9 



(38) 



where 



tan 29 = ^ < 9 < vr/2, (39) 

ill \ L2 

one can get the eigenstates for the one-exciton states je,) and the transition dipole moments 
D ei {i = 1, 2) corresponding to the transitions between the ground and single-exciton states 
as 
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u 



-1 



V a e2 J 



( \ 

A x cos 9 + A 2 sin 9 



\ 



-Ai sin 8 + A 2 cos 9 



(40) 



Here a £i = \ei),D e . and Ai = B^~\0),Di; Di and D 2 are the site transition moments. The 
two one-exciton energies are given by 



ft£t ei = HVt\ cos 9 + Ml 2 sin 2 9 + hJ sin 2$, 



hVt e2 = Kli sin 2 + KVL 2 cos 2 9 — hJ sin 2$, 



The two-exciton state wavefunction and its energy are as following 



(41) 



\e 3 ) = B+B+\0) = B+\0) 



(42) 



hQ e3 = hQi + Ml 2 



(43) 



The transition dipole moments between the single-exciton and two-exciton states are given 
by 



D eie3 = D\ sin 9 + D 2 cos 9, D e , zez = D\ cos 9 — D 2 sin 9 

However, the transition between the ground and two-exciton states is not allowed. 
In the eigenstate representation, the Hamiltonian of Eq. (l3"o"|) is rewritten as 

H= KQ ei -a ei )B+B ei -h J2 a eie .B+B ej + H bath - 



8=1,2,3 



+ B ei ) + IX ,„(!>>; IS,, + B+BJ] ■ E(t) 



i=l,2 



Here the interaction with the bath is given by 
/ \ 

h = u~ x H eb u = h 



\ 



^ a\ cos 2 9 + a 2 sin 2 9 \(ct 2 — «i) sin 29 ^ 



/ 



V 



\{a 2 — «i) sin 29 ot\ sin 2 9 + a 2 cos 2 # 
23 



/ 



(44) 



(45) 



(46) 



and 

a e3 — ot\ + a 2 , (47) 

for the single-exciton and two-exciton states, respectively. Eqs. (1461) and (14"T|) define the 
fluctuating parts of the single-exciton and two-exciton state transition frequencies. 

Consider various correlation functions. Assuming that baths acting on different chro- 
mophores are uncorrelated 

(a m (t)a n (0)) = 0ioim^n (48) 

and that the site energy fluctuation correlation functions are identical for the two monomers 
[551 155] , we get 



(a ei (*)<*«*! (0)) = (a e2 (f)a e2 (0)) = /T 2 K(t)(cos 4 9 + sin 4 6), 
(a e3 (t)a e3 (0))=2rt(t) (49) 

where K{t) = h~ 2 (ai(t)ai(0)) = h~ 2 (a2(t)a2(0)) = h~ 2 (a(t)a(0)). The further calculations 
simplify considerably if the off-diagonal part of the interaction with the bath in the exciton 
representation a eie2 = a e2ei in Eq. (l4T)i) can be neglected. This approximation is discussed in 
Refs.[55l[5g. 

The correlation function K(t) can be represented as the Fourier transform of the power 
spectrum $(cj) of hai(= ha.2) [59J 

/oo 
-00 

where 

$(-cj) = $(cj) exp(-Phu) (50) 
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Using Eq. (15171) . the real and imaginary parts of K(t) = K'(t) + iK"(t) can be written as 

poo 

K'{t) = / duj®(Lu)[l + exp(-(3huj)] cos cut 
Jo 

poo 

K"(t)= / du$(uj)[l-exp(-phu))smujt 
Jo 

In the high temperature limit one get 

POO 

K'(t) = 2 duj<$> (uj) cos cot 
Jo 

poo 

K"{t) = hf3 duj^(uj)ujsmujt 
Jo 

where K(0) = -^'(0) = 2 J °° dtu$(uj) = h 2 a 2 = hustP^ 1 ] °"2 an d cu st are a second cen- 
tral moment and the Stokes shift of the equilibrium absorption and luminescence spectra, 
respectively, for each monomer. 

Similar to SecEJ we will consider a = —u/h as a stochastic Gaussian variable with the 
correlation function (a(t)a(0)) = a 2 exp(— |£|/t s ) corresponding to the Gaussian-Markovian 
process. In this case the Fokker-Planck operators for the excited state of each monomer has 
the following form 



d 2 d 

+ {x-x m ) — + l) (51) 



\ ax z OX 

where x = quj\/j3 = a/^/a^ is a dimensionless generalized coordinate. Bearing in mind 
Eqs.f l49p . the Fokker-Planck operators for the eigenstates \j) = |0),|ej) of the exciton Hamil- 
tonian can be written by Eq.Q where x = 0, x ei = x e2 = x m (cos 4 #-|-sm 4 6) and x e3 = 2x m . 
The corresponding transition frequencies at the equilibrium nuclear coordinate of the ground 
electronic state are defined by Eqs. (l4Tj) and (1431) . 

Consider a homodimer complex consisting of identical molecules with Cli = Q 2 = fl and 
Di = D 2 = D. For this case, using Eqs.flU, (|10|1, (|4"I]h <P|) and (jiijh we obtain 9 = vr/4, 
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Figure 7: Populations of the ground (dotted line), single- (solid line) and two-exciton (dashed 
line) states of a homodimer complex after the completion of the pulse action as functions of 
$"(z/) for J = -300 cm- 1 (J < - J-aggregate), Q' = 2.9, t p0 = 10 fs, t s = 100 fs. The 
partial relaxation and the total models - lines with and without hollow circles, respectively. 



m ei2 = h(n ± j), m ea = 2hQ (52) 

D ei = D eie3 = V2D, D e2 = D e2e , = 

We thus need to consider only three states: |0), |ei) and |e 3 ), since state |e 2 ) is not excited 
with light. Letting |1), |2) and |3) represent |0), \e\) and |es), respectively, we arrive at 
a three-state system considered above where U21 = + J, u^i = 2fi, D 2 \ = D 32 = V%D, 

X\ 0, X 2 ~2^in) 2x m . 

FigH shows populations of single and two-exciton states after the excitation with a linear 
chirped pulse, Eqs.(!T5l) and ( fTBT) . as functions of $"(z/). Here the one-photon resonance for 
Franck-Condon transition 1 — » 2 occurs at the pulse maximum, i.e. u = U21 = Cl + J, and 
the Stokes shift of the equilibrium absorption and luminescence spectra for each monomer 
is equal to oj 1 ^ 011 = 400 cm -1 . Fig J7] also contrasts calculations using the total model (lines 
without hollow circles) with those of the partial relaxation model when only diagonal matrix 
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elements of the density matrix undergo diffusion (lines with hollow circles). FigJT] shows a 
good agreement between the calculation results for both models. 

Furthermore, one can see strong suppressing the population of the two-exciton state for 
negatively chirped (NC) pulse excitation. As a matter of fact, one can suppress or enhance 
two-exciton processes using positively or NC pulses. Our calculations (see table below) show 
twofold benefits of NC pulse excitation (<£>" = — 10 4 fs 2 ) with respect to the transform limited 
pulse ($" = 0) of the same duration (t p = 71 fs) and energy tuned to one-exciton transition: 
the population transfer to the single exciton state is larger, and that to the two-exciton state 
is smaller. 

Populations after the Transform limited 

NC pulse = -10 4 /s 2 , t p = 71 fs) 

completion of pulse action pulse (<&" = 0, t p — 71 fs) 

n 2 0.317 0.573 

n 3 0.208 0.057 

It is worthy to note good selective properties of chirped pulses, bearing in mind strong 

overlapping Franck-Condon transitions 1 — >• 2, lu 2 i, and 2^3, u 32 . Really, the corresponding 

frequencies differ by u 32 — ^21 = — 2J — |cj™ on for the model under consideration that 

comes to u 32 — ^21 = 300 cm -1 for the used values of parameters. On the other hand, 

the bandwidth of the absorption spectrum at half maximum for transition 2^3 comes to 

Au; = 2 a/2 In 2a|f « 1024 cmT x that is larger than u 32 — cv 2 i- Here off = (ft/?) -1 ^ 3 is 

the LF vibration contribution to a second central moment of an absorption spectrum for 

transition 2 — > 3 and u 2 ^ = {h[3)~ l {x 3 — x 2 ) 2 = \uj™ n is the corresponding Stokes shift. 

This issue can be understandable in terms of the competition between sequential and 

direct paths in a two-photon transition [36]. Consider a three-level atomic ladder system 

in the absence of relaxation with close transition frequencies uj 2 \ ~ uj 32 where uj 2 \ can be 

associated with one-exciton excitation and frequency co> 31 - with two-exciton excitation. The 
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system is affected by one phase modulated pulse of carrier frequency u, Eqs.(J2J), ffTBT) and 
( TT6|) . In the Appendix we have calculated the excited-state amplitude 03 due to two-photon 
transition 1 — 3 involving a nearly resonant intermediate level 2 for such system. Amplitude 
03 = cltp + < 1 s consists of two contributions. The first one a-pp corresponds exactly to that of 
the nonresonant two-photon transition. This contribution otp ~ and it is small 

for strongly chirped pulses [52] 

2|$»|»r p 2 (53) 

This result has a clear physical meaning. The point is that the phase structure (chirp) of the 
pulse determines the temporal ordering of its different frequency components. For a strongly 
chirped pulse when a pulse duration is much larger than that of the corresponding transform- 
limited one, one can ascribe to different instants of time the corresponding frequencies [52] . 
As a matter of fact, in the case under consideration different frequency components of the 
field are determined via values of the instantaneous pulse frequency u(t) for different instants 
of time. Therefore, only a small part of the whole pulse spectrum directly excites the two- 
photon resonance. 

The second contribution is given by [35] 

as = ~ D ^ h l llX E(lu 21 )E(lu 32 ) {1 - sgn[(u 21 - ^32)$" (w)]} (54) 

where E(Cj) is the Fourier transform of the positive frequency components of the field ampli- 
tude £ (t) exp[iipi (t)}. The consideration of the Appendix enables us to extend the results of 
Ref.[3E] to non-zero two-photon detuning Q 2 = ^31 — 2cj 7^ 0. Eq. (f51|) describes a sequential 
process, the contribution of which is a steplike function. This process can be suppressed 
when the pulse frequencies arrive in counter- intuitive order {u 32 before u 2 i) that occurs in 
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our simulations of a J-aggregate for NC excitation. FigJT] and the table above show that the 
selective properties of chirped pulses under discussion are conserved on strong field excita- 
tion and for broad transitions. The selective excitation of single and two-exciton states can 
be used for preparation of initial states for nonlinear spectroscopy based on pulse shaping 

pa eu. 

6 Strong interaction and STIRAP 

The three-state system under discussion enables us to consider STIRAP as well. STIRAP in 
molecules in solution was studied in Refs. [22] , where the solvent fluctuations were represented 
as a Gaussian random process, and in Ref . [30J , where the system-bath coupling was taken to 
be weak in the sense that the relaxation times were long in comparison to the bath correlation 
time, t c . Intense fields were shown in Ref. [30] to effectively slow down the dephasing when the 
energetic distance between the dressed (adiabatic) states exceeds l/r c . The point of the last 
paper is that in contrast to usual undressed states, which intersect, the dressed (adiabatic) 
states do not intersect. Therefore, the spectral density of the relaxation induced noise, 
which has a maximum at zero frequency, strongly diminishes for frequencies corresponding 
to the light-induced gap between dressed states, resulting in suppressing pure dephasing 
between the dressed states. In this section we show that this conclusion holds also for non- 
Markovian relaxation when the system-bath interaction is not weak and, therefore, can not 
be characterized only by r c . 

In the rotating wave approximation the Schrodinger equations for STIRAP in A-configuration 
can be written as follows 
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/ \ / 



dt 



a 2 
V° 3 ) 



\ ( \ 



I 



a 2 



U[ -Ml x /2 

-m x /2 u 2 -nn 2 /2 a 2 ( 55 ) 

I -hVL 2 /2 Ut 
where U' x = U x + fa^i and = U3 + hw 2 are "photonic replications" of effective parabolic 

potentials U x (x) and U 3 (x) (Eq. lflU]) ). respectively. We consider the two-photon resonance 

condition when U\ — uj 2 = (E3 — Ei)/h and x\ = x% = that would appear reasonable when 

|1) and |3) are different vibrational levels of the same electronic state. Then U[ = U^. 

Adiabatic states U ad corresponding to Eq. (l55l) can be found by equation 



/ 



det 



U[ - U ad -m x /2 

-hn x /2 U 2 - U ad -Ml 2 /2 

-m 2 /2 u'z - u ad J 

This gives the following adiabatic states 



V 



jjad _ Til _ Til 



u. 



ad 



~(u 2 + u[) ± -J(u 2 - u[) 2 + h 2 (nl + n 2 2 



(56) 



One can see that initial U[ and final U' z diabatic states coincide with one of adiabatis states 
Uq" 1 . For strong interaction the last will be well separated from other adiabatic states U± due 
to avoided crossing. Therefore, during STIRAP the system will remain in the same adiabatic 
state UQ d , which is U[ for t = —00 and for t = +00. Its evolution due to relaxation 
stimulated by LF vibrations can be described by the corresponding Fokker-Planck operator 
Lq^ = Li ;3 = Tg 1 (j^p + + 1 j describing diffusion in adiabatic potential UQ d = U[ = U^. 
This means that during transition 1 — > 3 the system motion along a generalized coordinate 
x does not change. In other words, such a transition will not be accompanied by pure 
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dephasing. This conclusion is a generalization of the previous result [30J relative to slowing 
down the dephasing in strong fields, which was obtained for weak system-bath interaction, 
to non-Markovian relaxation. 

7 Conclusion 

In this work we have studied the influence of ESA and two-exciton processes on a coherent 
population transfer with intense ultrashort chirped pulses in molecular systems in solution. 
An unified treatment of ARP in such systems has been developed using a three-state elec- 
tronic system with relaxation treated as a diffusion on electronic potential energy surfaces. 
We believe that such a simple model properly describes the main relaxation processes related 
to overdamped motions occurring in large molecules in solutions. 

Our calculations show that even with fast relaxation of a higher singlet state S n (n > 1) 
back to S±, ESA has a profound effect on coherent population transfer in complex molecules 
that necessitates a more accurate interpretation of the corresponding experimental data. In 
the absence of S n — > Si relaxation, the population of state |3), n 3 , strongly decreases when 
the chirp rate in the frequency domain |<3>"(z/)| increases. In order to appreciate the physical 
mechanism for such behavior, an approach to the total model - the relaxation-free model - 
was invoked. A comparison between the total model behavior and that of the relaxation-free 
model has shown that relaxation is responsible for strong decreasing n 3 as a function of 
$"(z/) in spite of meeting adiabatic criteria for both transitions 1 — ► 2 and 2 — > 3 separately. 
By this means usual criteria for ARP in a two-state system must be revised for a three-state 
system. 

To clarify this issue, we have developed a simple and physically clear model for ARP with 



31 



a linear chirped pulse in molecules with three electronic states in solution. The relaxation 
effects were considered in the framework of the LZ calculations putting in a third level 
generalized for random crossing of levels. The model has enabled us to obtain a simple 
formula for n^, Eq. (1331) . which is in excellent agreement with numerical calculations. In 
addition, the model gives us an extra criterion for coherent population transfer to those we 
have obtained before for a two-state system [27] . New criterion, Eq.( l34|) . implies conservation 
of the "counter-movement" of the "photonic repetitions" of states 1 and 3, in spite of random 
crossing of levels. 

Furthermore, we also applied our model to a molecular dimer consisting from two-level 
chromophores. A strong suppressing of two-exciton state population for NC pulse excitation 
of a J-aggregate has been demonstrated. We have shown that one can suppress or enhance 
two-exciton processes using positively or NC pulses. As a matter of fact, a method for 
quantum control of two-exciton states has been proposed. Our calculations show good 
selective properties of chirped pulses in spite of strong overlapping transitions related to the 
excitation of single- and two-exciton states. 

In the light of the limits [121 01] imposed on Eqs.([6]) and ([8]) for nondiagonal elements 
of the density matrix for the total model, we used a semiclassical (Lax) approximation 
(Eq-(JHD) (the partial relaxation model). The latter offers a particular advantage over the 
total model. The point is that the partial relaxation model can be derived not assuming 
the standard adiabatic elimination of the momentum for the non-diagonal density matrix, 
which is incorrect in the "slow modulation" limit [12]. A good agreement between calcu- 
lation results for the partial relaxation and the total models in the slow modulation limit 
(see Figsj2] and [7j) shows that a specific form of the relaxation term in the equations for 
nondiagonal elements of the density matrix p\2(x, t) and p23{x,t) is unimportant. By this 
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means the limits imposed on the last equation [13"1 HI] are of no practical importance for the 
problem under consideration in the slow modulation limit. This issue can be explained as 
follows. Our previous simulations [26J show that in spite of a quite different behavior of the 
coherences (nondiagonal density matrix elements) for the partial relaxation and the total 
models, their population wave packets pjj (x, t) behave much like. Since we are interested 
in the populations of the electronic states rij = J pjj (x, t) dx only, which are integrals of 
Pjj (x, t) over x, the distinctions between the two models under discussion become minimal. 

In conclusion, we have also demonstrated slowing down the pure dephasing on STIRAP 
in strong fields when the system-bath interaction is not weak (non-Markovian relaxation). 
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Appendix 

Consider a three-level system E\ < E 2 < E3 with close transition frequencies u 2 i ~ 
CU32 where uo 2 i can be associated with a single-exciton excitation and frequency uj 3 i - with 
two-exciton excitation. The system is affected by one phase modulated pulse of carrier 
frequency uj, Eq.([2]). The excited-state amplitude for a two-photon transition involving a 
nearly resonant intermediate level, can be written as 



D 32 D 21 
03 = — 7T 



E(u 21 )E(u 32 ) + -P dQ 
n J-00 



1 n f°° ^E(n + u)E{n 2 -q + lu) 



(57) 



n - (co>2i - Uj) 

where E(u) is the Fourier transform of the positive frequency components of the field am- 
plitude £ (t) exp[i(fi (t)}, ft = Co — uj, P is the principal Cauchy value, fl 2 = u 3 i — 2u is the 
two-photon detuning. For linear chirped excitation, Eqs. (115l) and (1161) . E(u) is given by 



E(Q) = ^or p0 exp{-ifi 2 [T p 2 /2 - i$"(u)]} (58) 
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Using Eq. (|58|) and introducing a new variable z = f2 — fi 2 /2, Eq. (j5"7|) can be written as 



Q3 = 



£>32£>2l7r 2 (£oTpo) 2 

2h 2 



{exp[--(5 2 + (fl a - 5) 2 )(r 2 /2 - 



1 ? Z" 00 

+ exp[--^(r p 2 /2-^"M)]-P/ 




} 



(59) 



where 5 = U21 — lu is one-photon detuning. The integral on the right-hand side of Eq. fl59|) can 
be evaluated for strongly chirped pulses [52], Eq. f[53l . when a pulse duration is much larger 
than that of the corresponding transform-limited one. In this case two frequency ranges give 
main contributions to the integral. The first one results from the method of stationary phase 
[53] . and it is localized near the two-photon resonance z = u — u> 3 i/2 = in the small range 
Aco> ~ 1/ In this case only a small part Au ~ 1/ y/\<&" (ojj \ of the whole pulse 
spectrum Au pu i se = 4:/r p0 directly excites the two-photon resonance, and the corresponding 
contribution ~ 1 / -y/|$"(u;)| is small due to Eq. (1531 . 

The second contribution to the integral is located near z = 5 — Q2/2 arid it is due to the 
pole at the real axes. This contribution is given by Eq.(|54l) of SecJHl 
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